Quantum correlation in three-qubit Heisenberg model with Dzyaloshinskii-Moriya interaction 
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We investigate the pairwise thermal quantum discord in a three-qubit XXZ model with Dzyaloshinskii-Moriya 
(DM) interaction. We find that the DM interaction can increase quantum discord to a fixed value in the anti- 
ferromagnetic system, but decreases quantum discord to a minimum first, then increases it to a fixed value in 
the ferromagnetic system. Abrupt change of quantum discord is observed, which indicates the abrupt change of 
groundstate. Dynamics of pairwise thermal quantum discord is also considered. We show that thermal discord 
vanishes in asymptotic limit regardless of its initial values, while thermal entanglement suddenly disappears at 
finite time. 
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I. INTRODUCTION 

Entanglement has been intensively investigated in recent 
years as it is thought to be a fundamental resource for quantum 
computational processing tasks, e.g., quantum computation, 
teleportation and dense coding . It is also considered to be 
the most nonclassical nature of quantum composite systems 
and impossible to be simulated within classical formalism [2]. 
Recently, it is found that there may be nonclassical correla- 
tions in the quantum states other than entanglement, which is 
fragile with respect to environment. Quantum discord (QD) 
is proposed to be a measure of quantumness of correlation 
of bipartite systems, which is arising from the difference be- 
tween two quantum extensions of the classical mutual infor- 
mation yfl. In the case of pure states, QD reduces to entropy 
of entanglement ||4J]. In mixed states, QD and entanglement 
are independent measures of correlation without simple rel- 
ative ordering [4, 5]. QD of disentangled state may not be 
zero 13]. Nonentanglement quantum discord may provide the 
quantum advantage in the deterministic quantum computation 
with one quantum bit (DQC1)|0, 0]. Though quantum dis- 
cord arises from quantum information, it may also be applied 
to condensed matter physics, such as indicating the quantum 
phase transition ifsl— Tlltl . 

Spin is thought to be a suitable candidate as qubits, and 
the Heisenberg model describes the basic spin-spin interac- 
tion lfl2ll . Naturally, this type of model has been extensively 
studied, and can be found in Refs. M. [I3TH6I1 . Some features 
of QD are revealed, and behavior of QD is characterized with 
tunable parameters, e.g., external field, temperature, spin-spin 
coupling constant etc. By changing the temperature and also 
by applying an external field, several remarkable effects for 
QD are observed, many of them in sharp contrast to the behav- 
ior observed for the entanglement. It is revealed that situations 
where QD increases with temperature while entanglement de- 
creases or while entanglement is zero in two-qubit Heisen- 
berg model QSQ . It is argued that models of low-dimensional 
magnetic materials may supplemented with Dzyaloshinsky- 
Moriya (DM) interaction, which is arising from the spin-orbit 



coupling fnll . The Dzyaloshinsky-Moriya exchange interac- 
tion describes superexchange between the interacting s pins , 
and is known to generate many dramatical features IT18l - l22ll . 
It is found that DM interaction can excite entanglement and 
teleportation fidelity. 

Decoherence is the main obstacle for practical implemen- 
tation of quantum computing. Dynamics of quantum discord 
in presence interaction with the environment can be found in 
Refs. 11231 - 12811 . Both Markovian and non-Markovian dynam- 
ics of QD were investigated theoretically and experimentally. 
It was revealed that quantum discord is more robust than en- 
tanglement against decoherence. In this sense, it is argued 
that the quantum algorithms based only on quantum discord 
is expected to be more robust than based on entanglement. Re- 
cently, Quantum correlation has been investigated experimen- 
tally 12911 . the sudden transition from a classical to a quantum 
decoherence regime is observed during the dynamics of a Bell 
diagonal state in a non-Markovian dephasing environment. 

In this paper we consider pairwise QD of three-qubit XXZ 
model with Dzyaloshinskii-Moriya interaction in a heat reser- 
voir. We start in Sec. [EI] for introducing three qubits XXZ 
model with DM interaction in a canonical ensemble. In 
Sec.|III]we investigate behavior of pairwise thermal QD with 
variation of system parameters, and present the respective nu- 
merical results. In Sec. |IV] We calculate the dissipative dy- 
namics of pairwise thermal quantum discord under Markovian 
environments. We analyze two types of quantum channel such 
as dephasing, depolarizing, by assuming each qubit coupling 
to local noisy independently. The conclusions are obtained in 
the last section. 



II. MODEL AND ITS SOLUTION 

The Hamiltonian for a XXZ model with z-component DM 
interaction can be expressed as 
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where cr%(a — x,y,z) denote the Pauli matrices for qubit n, J 
is the coupling constant, D is the z-component DM interaction 
strength, A is an anisotropy parameter, B is an external mag- 
netic field in the z-direction. Periodic boundary conditions are 
assumed as cr" = <t£. 



It is easy to calculate the eigenstates of the Hamiltonian: 



<h) = 1000), 

4>i) = 4=(|011> + |101> + |110», 
V3 

cf> 2 ) = -U|001> + I010) + |100», 
V3 



where 



4>4> 
4>s) 
4>6> 



-!-=[(* + V3)|011> + (i- V3)|101>-2flll0>], 
2V3 

— =[(i- V3)|011) + (;+ V3)|101)-2i|110)], 
2V3 

— =[(i+ V3)|001) + (j- V3)|010)-2/|100)], 
2V3 

— -[(i- V3)|001> + (i+ V3)|010)-2/|100)], 
2V3 

HID, (2) 



and 
27- 

£4 = 
E b -- 



the corresponding eigenvalues: £0 = jjaa -1- jd, c\ 



3 JA + 3B, E x = 
\JA-B,E 2 = 2J-\jA + B,E 3 = -B-J-\jA- V37D, 
- -B-J - \JA + V37D, E 5 = B-J - \ JA - V37D, 
B - J - \ JA + V37Z), £ 7 = § /A - 35, respectively. 



For a system in thermal equilibrium, the density matrix 
p(T) = exp(-BH)/Z, where B = l/ksT (£g is Boltzmann's 
constant, which we set it equal to 1) and Z = T r[exp(-BH)] 
is the partition function. Hence, in this model, we have 



y 

Z' 



3 exp (-3BB) + exp [0(2 J A - B)][exp (-2)87) 

+2 cosh( V3/37D) exp 037)] , 

3 exp (3BB) + exp &S(27A + B)][exp (-2/37) 

+2 cosh( V3/37Z)) exp 037)] , 

2 cosh(y6B) exp (2/37A)[exp (-2)87) 

+2 cosh( V3/37D) exp 087)], 

2 cosh03/i) exp (2/37A)[exp (-2)87) - 

cosh( V3/37D) exp 037) - / V3 sinh( yfepjD) exp 037)] , 

cosh(3y6fi) + cosh03B) exp (2/37 A) [exp (-2/37) 

+2 cosh( V3y67D) exp 037)], (5) 



and y* represents the complex conjugation of y. Quantum dis- 
cord of this state is called pairwise thermal quantum discord. 
In the paper we normalize the coupling constant between spins 
to 7 = -1 for the ferromagnetic case and 7=1 for the antifer- 
romagnetic case and study QD for the two cases separately. 



III. THERMAL QUANTUM DISCORD 

The quantum discord, which measures the quantum corre- 
lation of bipartite system pab, is defined as difference of two 
versions of quantum mutual information dU, 



where 



Q(Pab) = I(Pab) - C(p AB ), 



I(Pab) = S(p A ) + S(p B ) -S(pab), 



(6) 



(7) 



is the quantum mutual information, which measures the total 
correlation of system pab- S(p) = -Tr(plogp) is von Neu- 
mann entropy. Pa{B) is the reduced state of pab- And 



C(pab) = max{S(p A ) - / \ Pi$ (Pa)1* 



(8) 



P(T) = -YjZxpi-pEiMXM, 



(3) 



with Z = zZj^o exp(— j82?f). The two-qubit reduced density ma- 
trix is obtained by tracing all but the first two spins, namely, 
Pn(T) = 7>3[p(7)]. Due to the periodic boundary conditions, 
it is easy to check that all of the reduced matrices are equal. 
In the standard basis {|00>, |01), |10), |11)} the density matrix 
pu(T) is given by 



is the quantum conditional entropy, which is the maximum 
information of subsystem A obtained by performing a mea- 
surement on subsystem B measures classical correlation rt30ll . 
{W B } represents a set of von Neumann measurement on B. 
p' A = Tr B (n i B p AB W B )/Tr AB (n , B p AB 'n B ) is the state of A after 
obtaining outcome i on B, and p, = T rA B (W B pA B W B ). 

To compute QD between spins, we choose the set of pro- 
jectors {U' B = V\i)(i\V ] '\i = 0, 1, V e SU(2)} for subsystem B. 
Following the procedure of Ref. yD, we have the QD for the 
reduced thermal state as follow, 



p n (T) = 

H ' 6Z' 



u 0^ 

w y 

y* w 

v 



(4) 



QipniT)) = S( Pl ) + J^A k log 2 A k + min{S U S 2 }, (9) 



where p\ is the reduced state of p\ 2 (T), /U is the eigenvalue 
of the p l2 (T), Si = poKOo) + pih(6i) and S 2 = h(0 2 ) with 
h(ff) = — i=2 log 2 ^ ~ 1°S2 ^T' the parameters can be ex- 
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FIG. 1 : (Color online) thermal quantum discord as a function of D at 
the temperature T = 0.9, B = 0, (a) J = -1 and (b) J = 1 



FIG. 2: (Color online) thermal quantum discord as a function of A at 
the temperature T = 0.6, B = 0, (a) J = -1 and (b) J = 1 



pressed in terms of matrix elements as 6q = \(u - w)/(u + w)\, 
6>i = |(w - v)/(w + v)|, 6 2 = V("-v) 2 +4M 2 /6Z', p = 
(u + w)/6Z', and pi = (w + v)/6Z'. 

It is easy to verify that invariance of QD under D <-> -D, 
thus we just discuss the positive values of D. In Fig.Q] we plot 
QD as a function of DM interaction while other parameters are 
fixed (T = 0.9, B = 0). In the ferromagnetic case (J = — 1), as 
increasing D, QD decreases to a minimum value first and then 
increases to 0.3984. When D — > +oo, the density matrix © is 
given by p = j(\<f>4)(<f>4\ + \4>b){4>b\) for any finite temperature, 
the pairwise QD of this state is 0.3984. We notice that QD 
may indicate abrupt change of groundstate. When we set B = 
0, and A = 2, the groundstate depends on D, it may be 



r i(l0oX0ol + \<h)<<h\) 

|(|0O><0O| + I04><<N + 106X^,1 
<j(|04><*4l + \</>6)(<t>6\) 



\<h)<<h\) d = 



V3' 



De(^,+oo) 



(10) 



QD of the groundstate abruptly 
^- • n - 0.3333 ^ 



As tuning D from to 6 
changes around the critical point D — -4= : 
0.3984. Raising temperature mixes the groundstate with ex- 
cited state. Thus QD abruptly changes in the plot nearby the 
point D — 5/V3 « 2.8868. In the antiferromagnetic case 
(/ = 1), QD increases rapidly as D increasing, then tends to 
a fixed value. It is easy to check that the density matrix (f3j is 
5(103X031 + 105X051) at an Y finite temperature when D — > +oo, 
and the pairwise QD of this state is also equal to 0.3984. We 
observe that QD abruptly changes around two points. One can 
easily obtain the critical point of dashline, D — V3 and solid- 
line, D = where groundstate abruptly changes. Due to the 
temperature, they are not the exactly equal to the ones in the 
Fig. Eh). 

We now move to investigate how QD behaves as we change 



the anisotropy parameter at finite temperature. In Fig. 
we plot how QD depends on anisotropy parameter for T 
0.6 and zero magnetic field. Pairwise thermal entangle- 
ment of three-qubit Heisenberg model has been considered 
in Ref. l3~Ill . The authors found that thermal entanglement 
is absent from both the antiferromagnetic and the ferromag- 
netic XXZ model with anisotropy parameter A > 1 . However, 
there do exist nonzero values of QD in the vicinity of A > 1 . 
Once again we observe the abrupt change points in the plots. 
This behavior has been observed in iV-qubit XXZ model d 1 Qfl . 
It is argued that abrupt change of QD in the critical point of 
iV-qubit XXZ model is related to the maximization process 
of QD. The calculation shows that the groundstate abruptly 
changes at critical point A = 1 when D < y3. Hence we 
notice QD abruptly changes on the point A = 1 for D < y3. 
When D = 2, the critical point isA = V3-0.5~l.2321 near 
the abrupt change point of dot-dashline. The critical points 
we obtained may be not the exact point in the plot. This is 
because the energy gap between groundstate and first excited 
state is small. As depicted in Fig. |3a), the abrupt change of 
dotted line (D = 1.5) shift from the critical point, this is be- 
cause first excited state is close to the groundstate, the energy 
gap is AE = £4 - £0 ~ 0.4019. Comparing with ferromag- 
netic case, the QD of antiferromagnetic case behaves in a less 
different way as tuning the anisotropy parameter other than 
the abrupt change points are located on the different points. 
One can readily check that the abrupt change of QD is related 
to abrupt change of groundstate. Fig.|2]shows QD increases to 
fixed values with J A in both cases. The state of system would 
be |(|0o)(0o + 107X071) when JA — > -00, thus QD vanishes. 

We now turn to characterize the dependence of QD on mag- 
netic field B. It is easy to check that invariance of QD under 
magnetic filed B «-> —B. We just discuss the positive magnetic 
field B. Fig. [3fa)-3(b) show the behavior of QD of systems 
without DM interaction as changing B. The QD of ferromag- 
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FIG. 3: (Color online) thermal quantum discord as a function of B at 
temperature T = 0.9, (a) J = -1, D = 0, and (b) J = 1, D = 0, and 
(c)/=-l,D = 3, and(d)/= 1,£> = 3. 



FIG. 4: Dissipative dynamics of QD as a function of y at the temper- 
ature T = 0.5, B = 0, A = 1, (a) J = - 1 and (b) J = 1 for dephasing 
channel, and (c) J = -1, and (d) J = 1 for depolarizing channel. 



netic system may increase to a maximum, and then decrease. 
It is worth noting that magnetic field plays a beneficial role 
when JA > in the ferromagnetic system. However, in the an- 
tiferromagnetic system, magnetic field just slightly enhances 
QD when A = 0.5. The behavior of QD is plotted in Fig. [2c)- 
3(d) when DM interaction presents. By applying magnetic 
field, the value of QD supplemented with DM interaction may 
have a more larger range than the cases without DM interac- 
tion. However, QD finally vanishes as B increases, namely, 
strong magnetic field is harmful to QD. This is because den- 
sity matrix © is \(pi)((f>i\ (separable state) in this limit. 



IV. THERMAL STATE UNDER LOCAL DECOHERENCE 

Now we will focus on decoherence of thermal state when 
each qubit couples to local quantum noisy independently. The 
dynamics of three qubits interacting independently with indi- 
vidual environments is described by the solutions of the ap- 
propriate Born-Markov-Lindblad equations [32|, which can 
also be obtained by so called the Kraus operator approach [ 1 ] . 
Given the initial state p(0), its evolution equation can be writ- 
ten as 



by 



p(t) = Yj E ^.tP(W e 'Ut> 



(11) 



where the Kraus operators, E a jj T = E a ®Ep®E T |Q]] satisfy the 
completeness condition Yj a p T E \p T E afi.,T = I- The operators 
{E a } describe the quantum channel effects of one qubit. 

The dephasing channel, which describes the loss of quan- 
tum information without any energy dissipation QJJ] , is given 



E = 
Ei = 



1 

o 



o Vr 



(12) 



where y = 1 - e~ x ' and x denotes the decay rate associating 
with real physical process. To study pairwise QD, we trace out 
the 3rd site. The reduced density matrix of any pair of sites 
are equal because of the symmetry of the Hamiltonian and 
assumption of qubits coupling to environments independently. 
Taking the Eq. (0 as the initial state and using Eqs. (flTt and 
(ITZb . we have time-dependent state 



(13) 



where y and other parameters were defined above. We investi- 
gate the behavior of QD and entanglement in this state without 
magnetic filed as follow. Substituting Eq. (13) into Eq. (9), we 
can obtain the QD HH, 

M() + WQ M() + vv Uq Uq 

OZ/ o OZ/ o 





f " 





) 


1 





w 


(l-y)y 


6Z 7 





d-r)/ 


w 




lo 





v , 



s 



w + Z(l - y)\y \ w + 1(1 - y)\yo\ 
— log 2 - 



6Z; 



hi 



z=±i 

1 + W I +10 



6Z' 



(14) 



where uq, wq, Vq, yo and Z' represent the parameters defined in 
Eq. (0 for B — respectively, and 9 = max{|Mo - wq\/3ZL (1 - 
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FIG. 5: Dissipative dynamics of QD and concurrence as a function 
of y and D, at the temperature T = 0.5, B = 0, A = 0.5. (a) and (c) 
for QD; (b) and (d) for concurrence, (a) and (b) J = -1; (c) and (d) 
J = 1 




FIG. 6: Dissipative dynamics of QD and concurrence as a function 
of y and D, at the temperature T = 0.5, B = 0, A = 0.5. (a) and (c) 
for QD; (b) and (d) for concurrence, (a) and (b) J = -1; (c) and (d) 
J = 1 



y)\yQ\l3Z' }. To compare the behavior of entanglement with 
QD, we investigate the pairwise entanglement by using con- 
currence as the quantifier 113 311 . One can write out the concur- 
rence (C)forEq. (T3) JH, 

C = — max{0, (1 - y)\y \ - V^l- (15) 



Figure. |4ja)-4(b) show the numerical result for qubits cou- 
pling to dephasing channel. We notice that the derivative, 
dQ/dy, of dot-dash line (D = 3) in Fig.|4fa) abruptly changes 
around a point. The similar behavior is observed in the antif er- 
romagnetic system with nonzero D. The QD regardless of its 
initial value only disappear in the limit y — > 1. Fig.|3a)-5(d) 
shows the behavior of quantum discord and entanglement with 
variation of D and y for fixed temperature T - 0.5 and 
A = 0.5. QD decreases to zero when y = 1 (t — > oo) in 
both cases of J. Obviously, entanglement (C) decreases lin- 
early with y for fixed D in the plot Fig.|3b)-5(d), it undergoes 
a sudden death regardless of J and initial value of C at finite 
time. 



The depolarizing channel is an important type of quantum 
noise, it describes the process in which the density matrix is 
dynamically replaced by the state 7/2, 1 denoting identity ma- 
trix of a qubit. The Kraus operators simulate the effect of the 



depolarizing channel is given by [ 1] 




where y was defined above. We use the same procedure 
as above to evaluate QD and entanglement numerically. As 
shown in Fig. |5c)-4(d), QD decreases more fast under this 
type of noisy than dephasing channel. The derivative, dQ/dy, 
is continuous when qubits couple to this quantum channel, 
which is different from the case that qubit couples to dephas- 
ing channel. In Fig. HJa)-6(d) we plot the quantum discord 
and entanglement for the variation of D and y. When sys- 
tem couples to this kind of quantum noisy, entanglement once 
again disappears at finite time but QD does not. 



V. CONCLUSIONS 

In summary, We investigated the properties of thermal 
quantum discord within three-qubit Heisenberg model sup- 
plemented with DM interaction. We found that the DM in- 
teraction can both decrease and increase quantum correlation 
in the ferromagnetic XXZ system, but just increase quantum 
correlation rapidly in the antiferromagnetic XXZ system. The 
abrupt change of QD was observed as tuning parameters such 



5 



as D and A, which implies QD may signal the abrupt change 
of groundstate. The effects of magnetic fields on quantum dis- 
cord were also considered. Our result shows strong magnetic 
suppresses quantum correlation, while weak magnetic can in- 
crease or decrease QD by controlling J A. We utilized thermal 
state as initial condition to calculate the dynamics of pairwise 
QD under Markovian environments. We found pairwise en- 
tanglement of thermal state may occur sudden death but QD 
is robust enough to disappear in asymptotic time. 
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